Magnetic structure of azurite above the 1/3 magnetization plateau 
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The transition from tiie 1/3 magnetization plateau towards tlie saturation magnetization in azurite 
has been studied by low-temperature, high-magnetic-field, high-frequency proton nuclear magnetic 
resonance (NMR). The observed symmetrical splitting of the NMR spectra is incompatible with the 
longitudinal incommensurate order appearing when the longitudinal correlation function becomes 
dominant over the transverse one, which is the expected framework for the existence of the 2/3 
magnetization plateau. The spectra are rather interpreted in terms of a more standard transverse 
antiferromagnetic (canted) order. 
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Since the discovery of the magnetization plateau at 
1/3 of the saturation magnetization in the natural min- 
eral azurite [H, Cus (003)2(011)2, the system has been 
intensely studied as a model for the frustrated antifer- 
romagnetic Heisenberg spin-1/2 chain of "distorted di- 
amond" geometry Its large 1/3 plateau is of pure 
quantum origin, as it is based on the spin configuration 
where two more strongly coupled "dimer" spins are in 
a singlet state, while the third "monomer" spin is fully 
polarized 0, In this letter we focus on another in- 
teresting subject, the possible existence of a narrow 2/3 
plateau (so far unresolved by the magnetization mea- 
surements [11, 3) ill the middle of the increase of mag- 
netization from the 1/3 plateau towards the fully po- 
larized system 0, Q • Unlike the 1 /3 plateau, the 2 /3 
plateau is expected from Oshikawa's criterion fi\ to break 
the translational symmetry. Within the one-dimensional 
(ID) description (i.e., the Tomonaga-Luttinger descrip- 
tion), where the spin-spin correlation functions along the 
chain exhibit a power-law behavior, it is associated to 
a rather exceptional situation in which the longitudinal 
spin-spin correlation function becomes dominant over the 
transverse one due to the so-called "r? inversion" of the 
corresponding exponents (77 H < 77-'-) p, In this case 
an incommensurate (IC) longitudinal (i.e., parallel to the 
applied field) three-dimensional (3D) order is expected to 
be stabilized at low temperature, which can generate a 
plateau at the commensurate point. Theoretically it has 
been shown that the 2/3 plateau may exist in a diamond 
spin chain for certain values of exchange couplings 

Ji, J2, J3 lying within the relatively broad range of val- 
ues extracted from the experiments [H, Isj-fis!]. Namely, 
there is a controversy on these values, which seems to 
have converged only very recently [isj . 

Another, more conventional and thus less interesting 
situation appears when the transverse correlations are 
dominant, leading to a transverse antiferromagnetic (AF) 
order at low temperature, meaning that the dimer spins 



develop a canted polarization. The order would typically 
be of the Neel type where the direction of the transverse 
polarization is fixed, but it may also be a spiral IC or- 
der where the transverse polarization rotates along the 
chain. In this work we present nuclear magnetic res- 
onance (NMR) data which are clearly in favor of the 
Neel canted order. We thus exclude the more appeal- 
ing possibility of the longitudinal IC order leading to a 
2/3 plateau. 

In order to determine the microscopic nature of the 
spin structure of azurite above the 1/3 plateau, we per- 
formed a high-frequency (1.3-1.45 GHz) proton NMR at 
very high field (31-34 T) in the M9 magnet at LNCMI- 
Grenoble, using specially designed ("bottom tuned") 
NMR probe for the narrow bore •^He refrigerator. We 
mention that, in contrast to the previous NMR work at 
lower field [3j, the spin structure in this field range can- 
not be directly accessed by the copper NMR, because the 
signal is lost ("wiped out" by shortening of T2) due to en- 
hanced fluctuations of electronic spins. There is only one 
proton site in azurite, which however generates two NMR 
lines for an arbitrary direction of the magnetic field H, 
corresponding to two diamond chains per unit cell which 
are differently oriented with respect to H (see Fig. 1 in 
Ref. ^). When the field is applied within the ac-plane 
of mirror symmetry, both chains become equivalent and 
only one proton NMR line is observed. In this study the 
azurite single crystal was oriented with its c-axis approx- 
imately parallel to the magnetic field, so that the field is 
perpendicular to the chain 6-axis, and only one proton 
NMR line is expected for a homogeneous system. 

Indeed, in the field dependence of proton NMR spectra 
shown in Fig.[T]a single line is observed in the 1/3 plateau, 
that is below the second critical field 7?c2(0.6K) — 
31.35 T. The "activated" magnetic field dependence of 
the nuclear spin-lattice relaxation rate data {T^^ , not 
shown here) confirms that this field range corresponds to 
a linear closing of the gap at the end of the plateau, as is 
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FIG. 1: (color online) High-field proton NMR spectra of azu- 
rite covering the transition from the 1/3 plateau into the fully 
polarized system, taken at 0.6 K. Vertical offset of the spec- 
tra is proportional to the magnetic field, applied parallel to 
the c-axis. Frequency axis is taken relative to the proton 
Larmor frequency defined by its gyromagnetic ratio, ^7 = 
42.5774 MHz/T. 



also observed by high-field EPR [14|. At i7c2 the NMR 
line suddenly broadens and then splits into a nearly sym- 
metrical spectrum. Each side of the spectrum consists of 
a pair of lines just above i?c2- Each pair merges into a 
distribution with some structure, and transforms into a 
more regular Gaussian lineshape at higher fields. 

As the spectrum is nowhere completely split, it is dif- 
ficult to directly distinguish if the lineshape corresponds 
to a simple Neel type splitting into two (or 2x2) some- 
what overlapping lines centered at two frequencies /i and 
/2, or to an IC lineshape. In the latter case, for a ID si- 
nusoidal distribution of frequencies spreading between /i 
and /2, an NMR lineshape /(/) corresponds to the den- 
sity function Dif) cx [{f-fM - fT^^^ convoluted by 
some line-broadening function G{f), I{f) = D*G{f). In 
Fig.[5]we have performed the corresponding fits to two se- 
lected spectra presenting the most regular lineshape, one 
just above Hc2 and the other in the middle of the transi- 
tion towards the saturation. In the former case the best 
fit was achieved using ad hoc Lorentzian line-broadening 
function, while the latter spectra were better approxi- 
mated by a Gaussian line-broadening. In both cases the 



fits are in favor of an AF ordered system, while fits to 
the lineshape corresponding to an IC state clearly fail. 

Further quantitative analysis of the spectra is given in 
Fig. [21 We first observe that the first moment of the 
complete lineshape, which reflects an average spin polar- 
ization in the system, is barely changing although the 
magnetization of the system strongly increases. This can 
be attributed to a nearly zero value of the longitudinal 
coupling constant Azz of protons for this particular ori- 
entation of the sample, due to the cancellation of the hy- 
perfine (super-exchange induced) and the direct dipolar 
component. 

In order to quantify the splitting of the spectra shown 
in Fig. [1] which is apparently nearly field independent, 
we have fitted the spectra to two Gaussians (see Fig. [5]) 
and plotted in Fig. |3] their position taken with respect to 
the moment of each spectrum. For the three spectra 
close to Hc2, where the fit to 2x2 peaks is more appro- 
priate, the splitting is defined from the position of the 
weighted average of the two peaks on each side. In this 
way we find that the splitting of the spectra is highly 
symmetrical and nearly flat with very weakly expressed 
maximum at H2/3 — 32.6 T, close to the midpoint be- 
tween 1/3 magnetization and saturated magnetization. 
It turns out that the observed field dependence of the line 
splitting follows precisely the prediction (see Fig. ^ for 
the transverse staggered spin polarization 5^ in the sim- 
plest model of independent dimers, which has originally 
been employed to describe the Bose-Einstein condensa- 
tion (BEG) of triplets 15|, |l6[. We observe that close 



to Hc2 one of the Zeeman-split triplet states of a dimer 
becomes degenerate with the singlet. Mixing these two 
states produces a state where the two spins of a dimer ac- 
quire identical longitudinal polarizations S^- and opposite 
transverse polarizations 5'"'", which depend on the mixing 
coefficients ^5|| . The magnitudes of the two components 
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FIG. 2: (color online) Fits of two representative proton NMR 
spectra (solid lines) to the lineshapes corresponding to an AF 
ordered state (dotted lines) and to a ID IC state (dashed 
lines). 
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FIG. 3: (color online) Analysis of the spectra given in Fig. [T] 
The first moment, i.e. the average lineshift with respect to 
the Larmor frequency, in the 1/3 plateau (□), in the transition 
region (■), and close to full polarization (Kl). The position of 
two (•) or four (A, v) components of the nearly symmet- 
rically split spectra determined relative to the 1st moment. 
Close to the 1/3 plateau, where 2x2 components fit is ap- 
plied, • denotes the weighted average position of A and v- 
Solid and dashed lines are predicted field dependence for ID 
and 3D system of weakly interacting dimers, respectively (see 
the text). 



are related by an elliptical dependence: 



(1) 



where S'^ has a broad maximum at the point where the 
spins are half polarized, S^^^ = S-^iS^ = 1/4) = ^2/4, 
and develops very rapidly from zero, S^{S^ = 0, 1/2) = 
0, as soon as the system departs from zero or full polar- 
ization. In azurite the polarization of dimers in the 1 /3 
plateau is found to be 10% Q, so that we expect that 
at Hc2 the transverse polarization will develop starting 
right away from as much as S'^{S^ = 1/20) = 0.6 5*^^2.. 
Eq.[T]can be used to predict the field dependence of S*-*-, if 
(H) is known from the experiment or from the theory. 
High field magnetization data in azurite Q show that 
the increase of magnetization towards saturation is of the 
"arcsin" type corresponding to what is expected for a ID 
system. Our ID fit shown in Fig. [3] corresponds to the 
arcsin dependence defined by the initial point S^{Hc2) = 
1/20 and the midpoint 5*^(^^2/3) — 1/4. For compari- 
son, we also show a 3D fit, where the S'^{H) dependence 
is taken to be linear, passing through the same points. 
The ID fit is excellent, in particular when we take into 
account that the only fitting parameters are the over- 
all scaling defined by the observed maximum splitting 
of Afmax — 1.2 MHz and the corresponding field value 
i?2/3- The 3D fit is very close, and the experimental 
data do not really allow us to distinguish between the 
two fits. The fits confirm the description in terms of 
the transverse polarization, whose main characteristics 



is that strong increase of magnetization is accompanied 
by nearly field independent splitting of the spectra. We 
note that this particular type of behavior is not expected 
in a description in terms of an IC modulation of the lon- 
gitudinal polarization, which provides another argument 
against this latter interpretation. 

In the following we proceed by a quantitative analysis 
of the experimental value (A/mas) of the splitting. A 
transverse component of the spin polarization is shifting 
the NMR line through the off-diagonal (A^^;) elements of 
the coupling tensor. In general both the direct dipolar 
field, as well as the hyperfine field are expected to con- 
tribute to the proton coupling tensor. However, the lat- 
ter mechanism is based on an s-type wave function at the 
proton site, it is thus isotropic and cannot generate A^x 
elements. This means that the NMR splitting due to the 
transverse staggered S-^ is of pure dipolar origin and can 
be calculated for a given s pin order from the known crys- 
tallographic structure ITHlOf. In Fig. |4]we show so pre- 
dicted NMR linesplitting generated by ±S^^^ = ±\/2/4 
on dimer spins as a function of their direction in the plane 
perpendicular to the field applied along the c-axis. There 
are two possibilities for the stacking of S*^ (AFl and 
AF2) as regards the relative phase (0 or tt) of an AF spin 
arrangement in two chains in a unit cell. Comparing the 
prediction to the experimental value of ±1.2 MHz, we ob- 
serve that this is very close to the predicted ±1.35 MHz 
for the AFl type order with S'^ along the a* = a-axis. 
This particular direction agrees with the easy magnetiza- 
tion axis (perpendicular to the chains) , as defined by the 
low-field magnetization measurements [ij. However, the 
magnetization data correspond to the initial polarization 
of the monomer spins, while in the high-field measure- 
ment these spins are saturated and thus inert, and the 
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FIG. 4: (color online) Experimental value of the line split- 
ting (dashed line) compared to the predicted dipolar splitting 
for two possible AF orders (solid and dash-dotted line, as 
a function of the orientation of the ordered moments in the 
transverse plane. 
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active sites are dimers. Unfortunately, there is no direct 
information on the "focal" easy axis of dimers at high 
field. In general, azurite presents quite strong anisotropy 
of its magnetic properties, which has not yet been under- 
stood [11,1, SQ. 

One should remark that the predicted splitting from 
Eq. [T] corresponds to the maximum expected S"^, which 
is then to some extent reduced by thermal and quan- 
tum fluctuations. According to the complete mean field 
solution for the 3D interacting dimers 2l|, the ther- 



mal effects are negligible below 0.4 T™°^, where TJ"''^ is 
the maximum critical temperature. From measurements 
taken at 1.2 K we know that T™''^ is well above this 
temperature, meaning that the low temperature limit is 
fulfilled for our 0.6 K spectra. The effect of quantum 
fluctuations depends strongly on dimensionality and can 
be very important for purely ID systems. For exam- 
ple, in the CuBr4(C5Hi2N)2 (BPCB) spin ladder, where 
the effective 3D coupling is «50 times smaller than the 
inter-dimer interactions, the reduction of the transverse 
moment is ~50% 22, l23|. In azurite the 3D exchange 
couplings are estimated to be of the order of a few 
K, that is (only) one order of magnitude smaller than the 
(effective) ID couplings. At high magnetic fleld the ID 
Hamiltonian of azurite can be further reduced to the ef- 
fective XXZ, S =1/2 spin chain Hamiltonian describing 
only the dimer spins 0, H^. The corresponding cou- 
plings are of the order of 1 K, of the same size as the or- 
dering temperature, so that we can suspect that residual 
3D couplings are of comparable size. The 3D character 
of azurite at high field is thus probably important, mak- 
ing the effect of quantum fiuctuations relatively weak. 
For example, comparing our prediction for along the 
a*-axis to the experimental value we get the reduction 
of only salO%. We also recall that this description does 
not take into account the anisotropy terms of azurite's 
Hamiltonian [H, H, , which might also significantly 
modify the observed spin polarizations. 

Finally, within the effective XXZ Hamiltonian the 
anisotropy of the numerically optimized [25| effective cou- 
pling is found to be Jxy/ Jz =2.2, very close to the value 
of 2 obtained when the same description is a ppl ied to 
the spin ladder in the strong coupling regime 23]. For 
this latter system we know that in the 3D ordered phase 
the spins are canted and that no plateau is expected, in 
agreement with our NMR results. 

In conclusion, the high-field proton NMR spectra of 
azurite have been analyzed to determine which type of 
spin order is realized between the 1/3 magnetization 
plateau and the full polarization of the system. Both 
the observed line shape and the field dependence of the 
splitting of the spectra clearly correspond to the "conven- 
tional" transverse (canted) Neel order. They are both in- 
compatible with the longitudinal IC order, meaning that 
there is no 2/3 plateau in azurite, and that the corre- 
sponding "more interesting" physics does not apply to 



this system. These conclusions are supported by the re- 
cent theoretical estimates of the effective Hamiltonian for 
this system (isj. 
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